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With the advent of Fe-As based superconductivity it has become important to study how super¬ 
conductivity manifests itself in details of S7 Fe Mossbauer spectroscopy of conventional, Fe - bearing 
superconductors. To this end, the iron-based superconductor LugFesSis has been studied by 57 Fe 
Mossbauer spectroscopy over the temperature range from 4.4 K to room temperature with par¬ 
ticular attention to the region close to the superconducting transition temperature (T c = 6.1 K). 

Consistent with the two crystallographic sites for Fe in this structure, the observed spectra appear 
to have a pattern consisting of two doublets over the whole temperature range. The value of Debye 
temperature was estimated from temperature dependence of the isomer shift and the total spectral 
area and compared with the specific heat capacity data. Neither abnormal behavior of the hyperfine 
parameters at or near T c , nor phonon softening were observed. 

PACS numbers: 74.70.Dd, 76.80.+y 
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I. INTRODUCTION 

Since the discovery of Mossbauer effect spectroscopy 
more than half a century ago [IJ, superconductivity has 
been one of the states that has been investigated by 
this technique @]. Although Mossbauer spectroscopy is 
widely accepted as one of the most sensitive techniques in 
terms of energy resolution, it has not contributed signif¬ 
icant insight to studies of conventional superconductors. 
After the discovery of cuprate high temperature super¬ 
conductors (HTSC), Mossbauer spectroscopy was widely 
used for studies of these materials and reports of observa¬ 
tion of some anomalies in the spectral parameters in the 
vicinity of the superconducting critical temperature (T c ) 
[Ml were published. Due to the absence of commonly 
available Mossbauer nuclides in the cuprates, most stud¬ 
ies were accomplished either by partial substitution of 
copper atoms by 57 Fe and/or 119 Sn, or by using resonant 
isotopes of the rare earth metals, like 151 Eu, which in¬ 
creases the degree of difficulty of the measurements and 
reduces the clarity of the results sail. Recently, 
the discovery of iron-based superconductors, that natu¬ 
rally contain the common Mossbauer nuclide, 57 Fe, has 
triggered intense Mossbauer studies of these supercon¬ 
ductors [M3- Superconductivity in iron-pnictides is 
usually achieved by doping a magnetic parent compound 
with electrons or holes, or by application of chemical or 
physical pressure, and thereby, suppressing the magnetic 
order, suggesting that superconductivity and magnetism 
are closely related in this system. Although there are 
some studies on iron-based superconductors, which state 
that Mossbauer spectral parameters show anomalies near 
T c , in these materials it is hard to attribute the varia¬ 
tion of the hyperfine parameters observed by Mossbauer 
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spectroscopy to purely magnetic or purely superconduct¬ 
ing origins (Til .IfflIHi j . 

In order to study possible variations of hyperfine pa¬ 
rameters caused by the transition from the normal to the 
superconducting state in a conventional superconductor, 
we revisited lutetium-iron-silicide, Lu2Fe 3 Si5. Lu 2 Fe 3 Si5 
is a stoichiometric, Fe-containing, superconductor with 
relatively high T c « 6 K [l7i]. In a previous Mossbauer 
study [lj|, the non-magnetic nature of Fe was already 
confirmed. Hence, Lri2Fe 3 Si5 can be considered as an 
ideal compound to investigate the variation of hyperfine 
parameters caused only by the superconducting transi¬ 
tion without any complications associated by the absence 
of Mossbauer nucleus and/or the presence of magnetism. 
To the best of our knowledge no detailed, temperature 
dependent, 57 Fe Mossbauer spectroscopy measurements 
were performed on this material so far, and our goal is to 
shed some light on the applicability of Mossbauer spec¬ 
troscopy for studies of conventional, albeit multigap su¬ 
perconductors 0. 


II. EXPERIMENTAL DETAILS 

Polycrystalline samples of L^FeaSis were prepared by 
arc melting constituent elements with the nominal com¬ 
position of Lu2Fe3.32Si5.26 (corresponding to Lu 2 Fe 3 Si5 
+ Feo.32Sio.26) in Zr-gettered Ar atmosphere. Extra iron 
was added to suppress the formation of a Lu^FeSLt sec¬ 
ond phase and extra silicon was added to compensate 
apparent loss during the arc melting. To ensure the ho¬ 
mogeneity of the sample, the arc melting was repeated 
iteratively after flipping of the melted and resolidified in¬ 
got, for more than ten times. The weight loss was about 
0.26%. The arc-melted ingot was then sealed in an amor¬ 
phous silica tube, under a partial pressure of argon, and 
annealed at 1050 °C for 12 days. 
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Powder X-ray diffraction (XRD) was performed us¬ 
ing a Rigaku Miniflex diffractometer with Cu Ka radia¬ 
tion at room temperature (RT). The powder X-ray spec¬ 
tra of the samples were refined by Rietveld analysis us¬ 
ing the EXPGUI software [2C|. Magnetic measurements 
were performed using a Quantum Design Magnetic Prop¬ 
erty Measurement System SQUID magnetometer, spe¬ 
cific heat capacity was measured in a Quantum Design 
Physical Property Measurement System. 

Mossbauer spectroscopy measurements were per¬ 
formed using a SEE Co. conventional constant accelera¬ 
tion type spectrometer in transmission geometry with an 
57 Co(Rh) source, which had an initial intensity 50 mCi, 
kept at RT. The absorber was prepared in a powder form 
(10 mg of natural Fe/cm 2 ) by grinding of approximately 
175 mg piece of the arc-melted and annealed button. The 
absorber holder comprised two nested white Delrin cups. 
The powder was placed uniformly on the bottom of the 
larger cup and was held in place by a smaller cup. The 
absorber holder was locked in a thermal contact with a 
copper block with a temperature sensor and a heater, 
and aligned with the 7 - source and detector. The ab¬ 
sorber was cooled to a desired temperature using Janis 
model SHI-850-5 closed cycle refrigerator (with vibra¬ 
tions damping) that has long-term temperature stability 
better than 0.1 K at low temperature. The driver veloc¬ 
ity was calibrated by a-Fe foil and all isomer shifts (IS) 
are quoted relative to the a-Fe foil at RT. At first, spec¬ 
tra with maximum velocity 6 rnm/s and 3 rnrn/s were 
both measured at RT to check that no iron-containing 
impurity can be seen in the Mossbauer spectra. Then, 
three rounds of measurements, progressively focusing in 
on temperatures near T c = 6.1 K, were carried out: col¬ 
lecting 24 h with maximum velocity 2 nnn/s from 4.3 
K to 293.8 K (SI); collecting 48 h with maximum ve¬ 
locity 3 rnm/s from 4.4 K to 10 K (S2); collecting 48 
h with maximum velocity 3 rnm/s from 4.7 K to 6.4 K 
(S3). All the Mossbauer spectra were fitted by the com¬ 
mercial software package MossWinn [2lj |. in which the 
standard error of parameters can be estimated either by 
calculating and inverting the curvature matrix of the \ 2 
with respect to the fit parameters, or by the Monte Carlo 
method. The standard error of IS, quadroupole splitting 
(QS) and line width (T) in this work were obtained from 
the curvature matrix, while the error the area under the 
spectra here were obtained by Monte Carlo method by 
iterating 100 times. 


III. RESULTS AND DISCUSSION 
A. Structure and superconductivity 

The XRD pattern of the Ia^FeaSis polycrystalline 
sample is presented in Fig. 1. The majority of the peaks 
match to the tetragonal structure with the PA/mnc space 
group. The Rietveld refinement results in an estimate of 
~ 1.8(1) wt.% a-Fe impurity and a trace amount of un¬ 
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FIG. 1. (Color online) Rietveld refinement of the powder 
XRD spectrum of L^FesSis. Measured (black cross), calcu¬ 
lated intensities (red line) and background (green line) and 
difference curve (blue line) are shown. Vertical bars at the 
bottom indicate the positions of the Bragg reflections. The 
peaks of unknown phase are marked by the blue arrows. 
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FIG. 2. (Color online) Expanded view of the background part 
of the RT Lu 2 FesSi 5 Mossbauer spectrum in large velocity 
scale. Red arrows show the expected peaks positions for a- 
Fe.The inset is the full view of the spectrum. 


known phase. A 1.8 wt.% of a-Fe corresponds to 6.7 
% of Fe atoms in the a-Fe form in the ground sample, 
which, for un-enriched Fe, is below the resolution of 57 Fe 
Mossbauer spectroscopy to detect a-Fe. As can be seen 
in Fig. 2, there are no peaks associated with the a-Fe 
or with the other impurity in the Mossbauer spectrum, 
which means the unknown phase is either iron-free or, if 
contains iron, is below the resolution limit. Consequently, 
the Mossbauer spectra can be analyzed as a single - phase 
(Lu 2 Fe 3 Si 5 ) spectra. 

The superconductivity of the sample is confirmed by 
the dc susceptibility measurement in a magnetic field of 
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Temperature (K) 


10 Oe. As shown in Fig. 3 (a), the susceptibility data 
shows diamagnetic signal below ~ 6.2 K. The transition 
is sharp with a width of less than 0.4 K. Fig. 3 (b) shows 
the temperature dependent specific heat capacity (C p ). 
A sudden jump caused by superconducting transition can 
be observed below 6.3 K on cooling, the transition tem¬ 
perature is close to that obtained from the susceptibility 
data. Based on these two thermodynamic measurements 
we take T c = 6.1±0.1 K. 

Cp/T is also plotted as a function of T 2 in Fig. 4. A 
linear fit above the superconducting transition yields the 
values of y n (y n T is the electronic contribution to spe¬ 
cific heat capacity) and /3 n (/ 3 n T 3 represents the phonon 
contribution) of = 24.6(2) mJ/rnol K 2 , /3 n = 0.287(2) 
mJ/mol K 4 , which are very close to the previously re¬ 
ported values fzj] , From the /3 value, we can estimate the 
value of the Debye temperature (©d) using the relation: 
©d = (127r 4 A^rks/5/3) 1 / 3 , where N is Avogadro’s num¬ 
ber, r is the number of atoms per formula unit, and k b 
is the Boltzmann’s constant. We further obtained ©£> is 
408 K. The value of the normalized specific-heat jump at 
T c , AC/ 7 n T c , is ft: 1.06, a value that is smaller than the 
BCS value of 1.43, but consistent with the previously re¬ 
ported value 1.05 [22|. From the above characterizations, 
we can conclude that our sample is a bulk superconduc¬ 
tor and can be considered as single phase for Mossbauer 
measurements and analysis. 


B. Mossbauer results and discussion 


FIG. 3. (Color online) (a) Temperature dependent of zero 
field cooled magnetic susceptibility of Lu 2 Fe 3 Sis (H = 10 Oe); 
(b) the temperature dependent low temperature specific heat 
capacity. Vertical dashed line marks T c = 6.1 K. 



FIG. 4. (Color online) The specific heat capacity divided by 
temperature, Cp/T, as a function of T 2 for Lu 2 Fe 3 Sis. The 
dashed line represents the linear fit to the data in the normal 
state as described in the text. 


1. Symmetry of Fe sites and choice of the model 

In the Lu 2 Fe 3 Si 5 crystal structure, there are two, 
nonequivalent, Fe positions, Fei and Fen of 1 : 2 oc¬ 
cupation. Fei atoms are located at the 4d sites, which 
form ID chains along the c axis. Fen atoms are located 
at 8h sites, which form squares with planes perpendic¬ 
ular to the c axis. In each Fe position, the Fe atom is 
located in a polyhedron formed by Si atoms. Fei has four 
Si atoms at a distance of 2.31 A which form an irregular 
tetrahedron and two Si atoms at a distance of 2.54 A with 
the nearest Fe-Fe distance is 2.67 A. The Fen has four Si 
atoms at a distance of 2.34 A which are in the same face 
and form a quadrangle. On each side of the face, there 
are two Si atoms with the 2.35 A Fe-Si distances. The 
nearest Fe-Fe distance for Feu is 2.71 A. 

The anisotropic environments of the Fe atoms ensure 
nonzero electric field gradient (EFG) tensor at both sites. 
Hence, in analyzing the data, two doublets are expected. 
All the 57 Fe Mossbauer spectra of L^FeaSis over the 
whole observed temperature range share similar spectral 
shapes with a clear quadrupole splitting as shown in Fig. 
5 (a). A very small asymmetry and small shoulders were 
observed, which suggest that at least two subspectra are 
needed to resolve the spectra. However, due to the small 
splitting and consequently poor resolution of the spectra, 
the data can be analyzed with more than one set of pa- 












4 


TABLE I. Hyperfine parameters obtained by fitting using dif¬ 
ferent models as discussed in the text. IS is the isomer shift, 
QS is the quadrupole splitting and I is the relative area of the 
two subspectra. 


T 

K 

model 

site 

IS 

mm/s 

QS 

mm/s 

I 

% 

x 2 

RT 

1 

4d 

0.129(3) 

0.311(1) 

39.4 

1.63 



8h 

0.278(3) 

0.296(1) 

60.6 


RT 

2 

4d 

0.2157(9) 

0.513(6) 

36.4 

0.98 



8h 

0.2206(6) 

0.248(8) 

63.6 


4.4 

2 

4d 

0.320(2) 

0.523(4) 

52.8 

1.29 



8h 

0.330(1) 

0.202(3) 

47.2 



rameters. In the previous Mossbauer studies of R^FeaSis 
(R = rare earth), both (i) one doublet, and (ii) two dou¬ 
blets with fixed area ratio were employed to fit the spec¬ 
tra Hi]. In our approach to the fitting, two doublets 
without any restriction are used. In model 1, the two 
subspectra have close values of QS but obviously different 
values of IS, which is similar to the reported Sc 2 Fe 3 Sis fit 
result 0. In model 2, the two subspectra have similar 
values of IS, but distinct QS values. Two sets of param¬ 
eters yielding fits of acceptable quality can be obtained, 
the corresponding parameters of the RT spectrum fits are 
listed in the Table 1. The RT spectrum fitted using our 
two models is shown in Fig. 5 (b) and (c). The relative 
area of two subspectra in model 2 are closer to the the¬ 
oretical value 1 : 2 and \ 2 values are closer to 1, so we 
have chosen the model 2 to fit all the collected spectra. 
As an example, the fit using model 2 of the spectrum 
measured at 4.4 K from the S2 set is shown in Fig. 5 (d) 
and the corresponding parameters are listed in Table 1. 


2. Hyperfine parameters 

Fig. 6 summarizes the variation of hyperfine param¬ 
eters of the spectra with temperature. The plots in the 
right column present an expanded view of the low tem¬ 
perature range. 

The IS of Lu 2 Fe 3 Si 5 plotted as a function of tempera¬ 
ture is shown in Fig. 6 (a) and (b). The IS values of the 
4d and 8h sites at 294 K are 0.225(2) nun/s and 0.227(1) 
mm/s, respectively, which are slightly larger than the 
typical values for iron-silicon compounds in which Fe car¬ 
ries no moment [25j,[2(|, but are about 0.2 mm/s smaller 
than that in the iron-pnictide compounds [lOHT2i |. The 
temperature dependencies of the IS corresponding to the 
two sites are very similar, and no anomalies can be ob¬ 
served around T c (Fig. 6(b)). The IS values obtained 
from the fits includes contributions from both the chem¬ 
ical shift and the second-order Doppler shift, which is 
known to increase convexly upon decreasing temperature, 
due to gradual depopulation of the excited phonon states. 
However, it should be constant at low temperature, be¬ 
cause of the quantum mechanical zero-point motion. The 
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FIG. 5. (Color online) 57 Fe Mossbauer spectra of L^FesSis. 
(a)the spectrum at RT; the results of the fits using model 1 
and model 2 are shown in (b) and (c), respectively; (d) is the 
spectra collected at 4.4 K and fitted using the model 2. The 
model 1 and model 2 are described in detail in the text. 


chemical shift should not depend on temperature. The 
main contribution to this variation is from the second- 
order Doppler shift, which is usually described by Debye 
model: 


IS(T) = IS( 0) 


9 k B T T 3 /- 0d/t x 3 dx 

2Ah r W ) Jo e x — l' 


(1) 


where c is the velocity of light, M is the mass of the 57 Fe 
nucleus, and IS(0) is is the temperature-independent 
part, i.e. the chemical shift. A fit with Eq. (1) to the 
data of SI shown in Fig. 6 (a) yields 0£> = 517(18) K 
and 545(16) K for 4d and 8h sites, respectively. 

The quadrupole splittings in Lu 2 Fe 3 Si ,5 at the 4d and 
8h sites are 0.50(1) mm/s and 0.19(1) mm/s at 294 K, 
respectively. The magnitude of the QS is proportional 
to the 2 component V 22 of the EFG tensor, which is 
composed of two contributions: {V Z z)iig , from the lig- 
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FIG. 6. (Color online) Temperature dependence of parameters derived from fitting the of Lu 2 Fe 3 Sis Mossbauer spectra with 
model 2 as described in the text, (a) Isomer shift, (c) quadrupole splitting, (e) line width and (g) relative area of two subspectra. 
SI, S2, S3 - mark three measurement sets, the parameters are given for two Fe sites, 4d and 8h (see the text for more details). 
The lines in panel (a) are fits to the Eq. (1). The figures in the right column are the expansion one of the low temperature 
region of the corresponding left figure. The temperature stability for each of the measurements in the right column is better 
than that presented by the size of the symbols. 


and charges around the Mossbauer nucleus, and (V zz )val, 
from the valence electrons of Mossbauer nucleus. Usually, 
(y zz )ug is small and weakly dependent on temperature, 
whereas (V zz )vai is strongly temperature dependent. As 
can be seen from Fig. 5 (c). the QS of two sites are both 
almost temperature independent, which indicates the QS 
is mainly determined by the contribution from the ligand 


charge distribution around the o7 Fe sites. 

In order to get a better understanding of the elec¬ 
tronic origin of EFG at the Fe sites in L^FesSis, a 
first-principles calculation was performed using the full- 
potential linearized augmented plane wave method as 
embodied in the WIEN2K (27l. [28|. The generalized gra¬ 
dient approximation (GGA) suggested by Perdew, Burke, 
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Temperature (K) 


FIG. 7. (Color online) (a) Temperature dependence of the 
spectral area for SI set of measurements on a semi-log scale. 
The solid line is a fit to Eq.(2), as explained in the text, (b) 
Normalized to the values at « 6 K, temperature dependent 
spectral area data for all three sets of measurements. Inset: 
enlarged low temperature part. 


and Ernzerhof (PBE GGA) [29[ was employed for the 
exchange-correlation effects. Once the electron densities 
are calculated self-consistently and with high accuracy, 
the EFG tensor can be obtained from an integral over 
the non-spherical charge density. The principal compo¬ 
nent V 2 - for the 4d site is 2.25xl0 21 V/m 2 and asym¬ 
metry parameter r) = (V xx — V yy )/V zz = 0.457; the V zz 
for the 8h site is -1.16x 10 21 V/m 2 and r] = 0.307, which 
qualitatively agree with experimental results. It is found 
that the p-p and d-d interactions mainly contribute to 
the EFG of Lu 2 Fe 3 Si ,5 and p- electrons is play a signif¬ 
icant role for states far from the Fermi energy whereas 
the d-d interaction dominates around the Fermi energy, 
which is similar to what is found for the iron-pnictides 
superconductors mm 

The relative areas of the subspectra are determined by 
the proportion of the Fe atoms on different lattice sites. 
As mentioned above, for Lu 2 Fe 3 Si 5 , the theoretical rela¬ 
tive area of the two subspectra representing iron atoms 


at the 4d and 8h sites should be 33.3% and 66.7%. At 
RT the relative areas are close to expected value. How¬ 
ever, as shown in Fig. 6 (e) as temperature decreases, 
the relative area values deviate from theoretical values 
gradually, with the relative area for the 4d site becom¬ 
ing even larger than that for 8h site at low temperature. 
This phenomenon may be related to the variation of the 
linewidth, T, for the doublets corresponding to the two 
sites with reducing temperature. As shown in Fig. 6 (g), 
the T of 4d site increases slightly at lower temperature, 
but the T of 8h site is almost constant during the whole 
temperature range. The slightly increase of T of 4d site 
suggests the existence of instability at this lattice site at 
low temperature. 


3. The spectral area 

In Mossbauer studies of superconductivity, the varia¬ 
tion of the total spectral area has also been the focus 
of discussion. There are number of reports showing a 
decrease of the spectral area near T c due to the soften¬ 
ing of lattice with the opening of superconducting gap 
in cuprates but to the best of our knowledge only 
one report on iron-pnictide superconductors [32]. There 
are two kinds of behavior reported: a rapid decrease near 
T c 0,10; and a pit-like decrease either around T c or 
at higher temperature, serving as a precursor to T c Q. 
At the same time, some publications noted a poor repro¬ 
ducibility of those observations in the cuprates 0, [Til |. In 
our measurements, as can be seen in the Fig. 7, for the 
SI, we did not observe any abnormal variation around 
6.1 K. To make sure we didn’t miss any minor variation 
around T c , we remeasured the Mossbauer spectra of the 
same sample between 4.4 K to 10 K with higher density 
of the data points around 6.1 K for 48 h. Surprisingly, 
for data set S2, there is a sharp, 6%, decrease around 5.5 
K. Nevertheless, when we repeated the same measure¬ 
ment again, data set S3, this phenomenon disappeared. 
The sharp spectral area change seen in data set S2 is, in 
our opinion, most likely an artifact. The feature in data 
set S2 occurs resolvably below T c = 6.1 K with the very 
sharp and rather large jump occuring between spectra 
taken at 5.3 K and 5.7 K (with spectra at 5.7 K, 6.0 K, 
and 6.3 K all having normalized areas near 1.0). If this 
feature were associated with T 0 , it should have occured 
either between 5.7 K and 6.0 K or between 6.0 K and 
6.3 K. In addition, there are no corresponding anoma¬ 
lies in data set S2’s isomer shift, quadrupole splitting 
or relative areas (as shown in Fig. 6). A simple expla¬ 
nation for such behavior could be, among others, some 
mechanical shift/rearrangement of the powder compos¬ 
ing the absorber. This observation gives a warning that 
even in the measurement of iron-containing stoichiomet¬ 
ric material, irreproducibilities/artifacts might exist and 
should be addressed appropriately. 

Finally, we also fitted the temperature dependence of 
area under the two doublets line of SI measurements with 
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Debye model: 


/ 


-3 E 2 I T 2 /- eD/T 


xdx 
e x - 1 


}], ( 2 ) 


where / is the recoilless fraction, which is proportional 
to the area for thin sample and E 7 is the 7 -ray energy. 
This expression also allows to estimate the value of 0£>. 
We obtained the 0^ = 385(8) K which is very close to 
the 408 K obtained from the analysis of the low temper¬ 
ature specific heat capacity data and about 140 K less 
than the value estimated by temperature dependence of 
IS. A similar difference was found earlier in studies of e.g. 
FeSeQ.sTeo.s and 57 Fe- doped YBa 2 Cu 306.8 compounds 
[32l . I34| . This discrepancy may be explained by the fact 
the area reflects the average mean-square displacements, 
while IS related to the mean-square velocity of the Moss- 
bauer atom. Both quantities may respond in different 
ways to the lattice anharmonicities. 


principles calculations yield the values of EFG that qual¬ 
itatively agree with the experiment. The Debye temper¬ 
ature was estimated by the temperature dependence of 
specific heat capacity, spectral area and IS. The @25 ob¬ 
tained from temperature dependence of spectral area and 
heat capacity are very similar, but about 140 K smaller 
than the value estimated by the IS variation with tem¬ 
perature. Additionally, we didn’t observe any obvious, 
abnormal variation of hyperfme parameters around T c . 
Two possibilities could lead to this result: the opening 
of the superconducting gap doesn’t bring variation of the 
environment at Fe site at all; the T c of this system is too 
low, and Mossbauer spectroscopy is not sensitive enough 
to detect the minute change. 
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